Membrane specificity of the human cholesterol transfer protein STARD4.

STARD4 regulates cholesterol homeostasis by transferring cholesterol between the plasma membrane and endoplasmic reticulum. The STARD4 structure features a helix-grip fold surrounding a large hydrophobic cavity holding the sterol. Its access is controlled by a gate formed by the Ω1 and Ω4 loops and the C-terminal α-helix. Little is known about the mechanisms by which STARD4 binds to membranes and extracts/releases cholesterol. All available structures of STARD4 are without a bound sterol and display the same closed conformation of the gate. The cholesterol transfer activity of the mouse STARD4 is enhanced in the presence of anionic lipids, and in particular of phosphatidylinositol biphosphates (PIP2) for which two binding sites were proposed on the mouse STARD4 surface. Yet only one of these sites is conserved in human STARD4. We here report the results of a liposome microarray-based assay and microseconds-long molecular dynamics simulations of human STARD4 with complex lipid bilayers mimicking the composition of the donor and acceptor membranes. We show that the binding of apo form of human STARD4 is sensitive to the presence of PIP2 through two specific binding sites, one of which was not identified on mouse STARD4. We report two novel conformations of the gate in holo-STARD4: a yet-unobserved close conformation and an open conformation of Ω4 shedding light on the opening/closure mechanism needed for cholesterol uptake/release. Overall, the modulation of human STARD4 membrane-binding by lipid composition, and by the presence of the cargo supports the capacity of human STARD4 to achieve directed transfer between specific organelle membranes.

Overlay databank unlocks data-driven analyses of biomolecules for all.

Tools based on artificial intelligence (AI) are currently revolutionising many fields, yet their applications are often limited by the lack of suitable training data in programmatically accessible format. Here we propose an effective solution to make data scattered in various locations and formats accessible for data-driven and machine learning applications using the overlay databank format. To demonstrate the practical relevance of such approach, we present the NMRlipids Databank-a community-driven, open-for-all database featuring programmatic access to quality-evaluated atom-resolution molecular dynamics simulations of cellular membranes. Cellular membrane lipid composition is implicated in diseases and controls major biological functions, but membranes are difficult to study experimentally due to their intrinsic disorder and complex phase behaviour. While MD simulations have been useful in understanding membrane systems, they require significant computational resources and often suffer from inaccuracies in model parameters. Here, we demonstrate how programmable interface for flexible implementation of data-driven and machine learning applications, and rapid access to simulation data through a graphical user interface, unlock possibilities beyond current MD simulation and experimental studies to understand cellular membranes. The proposed overlay databank concept can be further applied to other biomolecules, as well as in other fields where similar barriers hinder the AI revolution.

Investigating Polypharmacology through Targeting Known Human Neutrophil Elastase Inhibitors to Proteinase 3.

Using a combination of multisite λ-dynamics (MSλD) together with in vitro IC50 assays, we evaluated the polypharmacological potential of a scaffold currently in clinical trials for inhibition of human neutrophil elastase (HNE), targeting cardiopulmonary disease, for efficacious inhibition of Proteinase 3 (PR3), a related neutrophil serine proteinase. The affinities we observe suggest that the dihydropyrimidinone scaffold can serve as a suitable starting point for the establishment of polypharmacologically targeting both enzymes and enhancing the potential for treatments addressing diseases like chronic obstructive pulmonary disease.

Dissecting peripheral protein-membrane interfaces.

Peripheral membrane proteins (PMPs) include a wide variety of proteins that have in common to bind transiently to the chemically complex interfacial region of membranes through their interfacial binding site (IBS). In contrast to protein-protein or protein-DNA/RNA interfaces, peripheral protein-membrane interfaces are poorly characterized. We collected a dataset of PMP domains representative of the variety of PMP functions: membrane-targeting domains (Annexin, C1, C2, discoidin C2, PH, PX), enzymes (PLA, PLC/D) and lipid-transfer proteins (START). The dataset contains 1328 experimental structures and 1194 AphaFold models. We mapped the amino acid composition and structural patterns of the IBS of each protein in this dataset, and evaluated which were more likely to be found at the IBS compared to the rest of the domains' accessible surface. In agreement with earlier work we find that about two thirds of the PMPs in the dataset have protruding hydrophobes (Leu, Ile, Phe, Tyr, Trp and Met) at their IBS. The three aromatic amino acids Trp, Tyr and Phe are a hallmark of PMPs IBS regardless of whether they protrude on loops or not. This is also the case for lysines but not arginines suggesting that, unlike for Arg-rich membrane-active peptides, the less membrane-disruptive lysine is preferred in PMPs. Another striking observation was the over-representation of glycines at the IBS of PMPs compared to the rest of their surface, possibly procuring IBS loops a much-needed flexibility to insert in-between membrane lipids. The analysis of the 9 superfamilies revealed amino acid distribution patterns in agreement with their known functions and membrane-binding mechanisms. Besides revealing novel amino acids patterns at protein-membrane interfaces, our work contributes a new PMP dataset and an analysis pipeline that can be further built upon for future studies of PMPs properties, or for developing PMPs prediction tools using for example, machine learning approaches.

Phosphatidylcholine Cation-Tyrosine π Complexes: Motifs for Membrane Binding by a Bacterial Phospholipase C.

Phosphatidylinositol-specific phospholipase C (PI-PLC) enzymes are a virulence factor in many Gram-positive organisms. The specific activity of the Bacillus thuringiensis PI-PLC is significantly increased by adding phosphatidylcholine (PC) to vesicles composed of the substrate phosphatidylinositol, in part because the inclusion of PC reduces the apparent Kd for the vesicle binding by as much as 1000-fold when comparing PC-rich vesicles to PI vesicles. This review summarizes (i) the experimental work that localized a site on BtPI-PLC where PC is bound as a PC choline cation-Tyr-π complex and (ii) the computational work (including all-atom molecular dynamics simulations) that refined the original complex and found a second persistent PC cation-Tyr-π complex. Both complexes are critical for vesicle binding. These results have led to a model for PC functioning as an allosteric effector of the enzyme by altering the protein dynamics and stabilizing an 'open' active site conformation.

Standard Binding Free Energy and Membrane Desorption Mechanism for a Phospholipase C.

Peripheral membrane proteins (PMPs) bind temporarily to cellular membranes and play important roles in signaling, lipid metabolism, and membrane trafficking. Obtaining accurate membrane-PMP affinities using experimental techniques is more challenging than for protein-ligand affinities in an aqueous solution. At the theoretical level, calculation of the standard protein-membrane binding free energy using molecular dynamics simulations remains a daunting challenge owing to the size of the biological objects at play, the slow lipid diffusion, and the large variation in configurational entropy that accompanies the binding process. To overcome these challenges, we used a computational framework relying on a series of potential-of-mean-force (PMF) calculations including a set of geometrical restraints on collective variables. This methodology allowed us to determine the standard binding free energy of a PMP to a phospholipid bilayer using an all-atom force field. Bacillus thuringiensis phosphatidylinositol-specific phospholipase C (BtPI-PLC) was chosen due to its importance as a virulence factor and owing to the host of experimental affinity data available. We computed a standard binding free energy of -8.2 ± 1.4 kcal/mol in reasonable agreement with the reported experimental values (-6.6 ± 0.2 kcal/mol). In light of the 2.3-µs separation PMF calculation, we investigated the mechanism whereby BtPI-PLC disengages from interactions with the lipid bilayer during separation. We describe how a short amphipathic helix engages in transitory interactions to ease the passage of its hydrophobes through the interfacial region upon desorption from the bilayer.

Specificity of Loxosceles α clade phospholipase D enzymes for choline-containing lipids: Role of a conserved aromatic cage.

Spider venom GDPD-like phospholipases D (SicTox) have been identified to be one of the major toxins in recluse spider venom. They are divided into two major clades: the α clade and the ß clade. Most α clade toxins present high activity against lipids with choline head groups such as sphingomyelin, while activities in ß clade toxins vary and include preference for substrates containing ethanolamine headgroups (Sicarius terrosus, St_ßIB1). A structural comparison of available structures of phospholipases D (PLDs) reveals a conserved aromatic cage in the α clade. To test the potential influence of the aromatic cage on membrane-lipid specificity we performed molecular dynamics (MD) simulations of the binding of several PLDs onto lipid bilayers containing choline headgroups; two SicTox from the α clade, Loxosceles intermedia αIA1 (Li_αIA) and Loxosceles laeta αIII1 (Ll_αIII1), and one from the ß clade, St_ßIB1. The simulation results reveal that the aromatic cage captures a choline-headgroup and suggest that the cage plays a major role in lipid specificity. We also simulated an engineered St_ßIB1, where we introduced the aromatic cage, and this led to binding with choline-containing lipids. Moreover, a multiple sequence alignment revealed the conservation of the aromatic cage among the α clade PLDs. Here, we confirmed that the i-face of α and ß clade PLDs is involved in their binding to choline and ethanolamine-containing bilayers, respectively. Furthermore, our results suggest a major role in choline lipid recognition of the aromatic cage of the α clade PLDs. The MD simulation results are supported by in vitro liposome binding assay experiments.

Phospholipids in Motion: High-Resolution 31P NMR Field Cycling Studies.

Diverse phospholipid motions are key to membrane function but can be quite difficult to untangle and quantify. High-resolution field cycling 31P NMR spin-lattice relaxometry, where the sample is excited at high field, shuttled in the magnet bore for low-field relaxation, then shuttled back to high field for readout of the residual magnetization, provides data on phospholipid dynamics and structure. This information is encoded in the field dependence of the 31P spin-lattice relaxation rate (R1). In the field range from 11.74 down to 0.003 T, three dipolar nuclear magnetic relaxation dispersions (NMRDs) and one due to 31P chemical shift anisotropy contribute to R1 of phospholipids. Extraction of correlation times and maximum relaxation amplitudes for these NMRDs provides (1) lateral diffusion constants for different phospholipids in the same bilayer, (2) estimates of how additives alter the motion of the phospholipid about its long axis, and (3) an average 31P-1H angle with respect to the bilayer normal, which reveals that polar headgroup motion is not restricted on a microsecond timescale. Relative motions within a phospholipid are also provided by comparing 31P NMRD profiles for specifically deuterated molecules as well as 13C and 1H field dependence profiles to that of 31P. Although this work has dealt exclusively with phospholipids in small unilamellar vesicles, these same NMRDs can be measured for phospholipids in micelles and nanodisks, making this technique useful for monitoring lipid behavior in a variety of structures and assessing how additives alter specific lipid motions.

Martini 3: a general purpose force field for coarse-grained molecular dynamics.

The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

Classification and phylogeny for the annotation of novel eukaryotic GNAT acetyltransferases.

The enzymes of the GCN5-related N-acetyltransferase (GNAT) superfamily count more than 870 000 members through all kingdoms of life and share the same structural fold. GNAT enzymes transfer an acyl moiety from acyl coenzyme A to a wide range of substrates including aminoglycosides, serotonin, glucosamine-6-phosphate, protein N-termini and lysine residues of histones and other proteins. The GNAT subtype of protein N-terminal acetyltransferases (NATs) alone targets a majority of all eukaryotic proteins stressing the omnipresence of the GNAT enzymes. Despite the highly conserved GNAT fold, sequence similarity is quite low between members of this superfamily even when substrates are similar. Furthermore, this superfamily is phylogenetically not well characterized. Thus functional annotation based on sequence similarity is unreliable and strongly hampered for thousands of GNAT members that remain biochemically uncharacterized. Here we used sequence similarity networks to map the sequence space and propose a new classification for eukaryotic GNAT acetyltransferases. Using the new classification, we built a phylogenetic tree, representing the entire GNAT acetyltransferase superfamily. Our results show that protein NATs have evolved more than once on the GNAT acetylation scaffold. We use our classification to predict the function of uncharacterized sequences and verify by in vitro protein assays that two fungal genes encode NAT enzymes targeting specific protein N-terminal sequences, showing that even slight changes on the GNAT fold can lead to change in substrate specificity. In addition to providing a new map of the relationship between eukaryotic acetyltransferases the classification proposed constitutes a tool to improve functional annotation of GNAT acetyltransferases.

Dynamics-function relationship in the catalytic domains of N-terminal acetyltransferases.

N-terminal acetyltransferases (NATs) belong to the superfamily of acetyltransferases. They are enzymes catalysing the transfer of an acetyl group from acetyl coenzyme A to the N-terminus of polypeptide chains. N-terminal acetylation is one of the most common protein modifications. To date, not much is known on the molecular basis for the exclusive substrate specificity of NATs. All NATs share a common fold called GNAT. A characteristic of NATs is the ß6ß7 hairpin loop covering the active site and forming with the α1α2 loop a narrow tunnel surrounding the catalytic site in which cofactor and polypeptide meet and exchange an acetyl group. We investigated the dynamics-function relationships of all available structures of NATs covering the three domains of Life. Using an elastic network model and normal mode analysis, we found a common dynamics pattern conserved through the GNAT fold; a rigid V-shaped groove formed by the ß4 and ß5 strands and splitting the fold in two dynamical subdomains. Loops α1α2, ß3ß4 and ß6ß7 all show clear displacements in the low frequency normal modes. We characterized the mobility of the loops and show that even limited conformational changes of the loops along the low-frequency modes are able to significantly change the size and shape of the ligand binding sites. Based on the fact that these movements are present in most low-frequency modes, and common to all NATs, we suggest that the α1α2 and ß6ß7 loops may regulate ligand uptake and the release of the acetylated polypeptide.

The Arabidopsis (ASHH2) CW domain binds monomethylated K4 of the histone H3 tail through conformational selection.

Chromatin post-translational modifications are thought to be important for epigenetic effects on gene expression. Methylation of histone N-terminal tail lysine residues constitutes one of many such modifications, executed by families of histone lysine methyltransferase (HKMTase). One such protein is ASHH2 from the flowering plant Arabidopsis thaliana, equipped with the interaction domain, CW, and the HKMTase domain, SET. The CW domain of ASHH2 is a selective binder of monomethylation at lysine 4 on histone H3 (H3K4me1) and likely helps the enzyme dock correctly onto chromatin sites. The study of CW and related interaction domains has so far been emphasizing lock-key models, missing important aspects of histone-tail CW interactions. We here present an analysis of the ASHH2 CW-H3K4me1 complex using NMR and molecular dynamics, as well as mutation and affinity studies of flexible coils. ß-augmentation and rearrangement of coils coincide with changes in the flexibility of the complex, in particular the η1, η3 and C-terminal coils, but also in the ß1 and ß2 strands and the C-terminal part of the ligand. Furthermore, we show that mutating residues with outlier dynamic behaviour affect the complex binding affinity despite these not being in direct contact with the ligand. Overall, the binding process is consistent with conformational selection. We propose that this binding mechanism presents an advantage when searching for the correct post-translational modification state among the highly modified and flexible histone tails, and also that the binding shifts the catalytic SET domain towards the nucleosome. DATABASES: Structural data are available in the PDB database under the accession code 6QXZ. Resonance assignments for CW42 in its apo- and holo-forms are available in the BMRB database under the accession code 27251.

Capturing Choline-Aromatics Cation-π Interactions in the MARTINI Force Field.

Cation-π interactions play an important role in biomolecular recognition, including interactions between membrane phosphatidylcholine lipids and aromatic amino acids of peripheral proteins. While molecular mechanics coarse grain (CG) force fields are particularly well suited to simulate membrane proteins in general, they are not parameterized to explicitly reproduce cation-π interactions. We here propose a modification of the polarizable MARTINI coarse grain (CG) model enabling it to model membrane binding events of peripheral proteins whose aromatic amino acid interactions with choline headgroups are crucial for their membrane binding. For this purpose, we first collected and curated a dataset of eight peripheral proteins from different families. We find that the MARTINI CG model expectedly underestimates aromatics-choline interactions and is unable to reproduce membrane binding of the peripheral proteins in our dataset. Adjustments of the relevant interactions in the polarizable MARTINI force field yield significant improvements in the observed binding events. The orientation of each membrane-bound protein is comparable to reference data from all-atom simulations and experimental binding data. We also use negative controls to ensure that choline-aromatics interactions are not overestimated. We finally check that membrane properties, transmembrane proteins, and membrane translocation potential of mean force (PMF) of aromatic amino acid side-chain analogues are not affected by the new parameter set. This new version "MARTINI 2.3P" is a significant improvement over its predecessors and is suitable for modeling membrane proteins including peripheral membrane binding of peptides and proteins.

Interfacial Aromatics Mediating Cation-π Interactions with Choline-Containing Lipids Can Contribute as Much to Peripheral Protein Affinity for Membranes as Aromatics Inserted below the Phosphates.

Membrane-binding interfaces of peripheral proteins are restricted to a small part of their exposed surface, so the ability to engage in strong selective interactions with membrane lipids at various depths in the interface, both below and above the phosphates, is an advantage. Driven by their hydrophobicity, aromatic amino acids preferentially partition into membrane interfaces, often below the phosphates, yet enthalpically favorable interactions with the lipid headgroups, above the phosphate plane, are likely to further stabilize high interfacial positions. Using free-energy perturbation, we calculate the energetic cost of alanine substitution for 11 interfacial aromatic amino acids from 3 peripheral proteins. We show that the involvement in cation-π interactions with the headgroups (i) increases the ΔΔGtransfer as compared with insertion at the same depth without cation-π stabilization and (ii) can contribute at least as much as deeper insertion below the phosphates, highlighting the multiple roles of aromatics in peripheral membrane protein affinity.

Peptidomimetic inhibitors targeting the membrane-binding site of the neutrophil proteinase 3.

Proteinase 3 (PR3), together with other serine proteases, such as neutrophil elastase (NE) and cathepsin G (CG), regulates inflammatory and immune responses. However, in comparison with NE and CG, there is increasing evidence that PR3 functions significantly differ. In particular, PR3 can bind to cell membranes and such membrane-bound PR3 (mbPR3) might be differently involved in the activation of cytokines, growth factors, cellular receptors, and in the regulation of cell apoptosis. For instance, PR3 membrane binding can block some "eat me" signals, notably, phosphatidylserine membrane lipid, and facilitate non-resolving inflammation. Based on the clear evidence that PR3 membrane binding affects the biological functions of PR3, we designed peptidomimetic inhibitors that can remove mbPR3 from the membrane surface in vitro without influencing PR3 catalytic activity. Such inhibitors, which specifically target PR3 binding to membranes, are still lacking. In particular, we found peptidomimetics that inhibit binding of PR3 to POPC:PS liposomes, which mimic the biological environment of PR3.

Cation-π Interactions between Methylated Ammonium Groups and Tryptophan in the CHARMM36 Additive Force Field.

Cation-π interactions between tryptophan and choline or trimethylated lysines are vital for many biological processes. The performance of the additive CHARMM36 force field against target quantum mechanical data is shown to reproduce QM equilibrium geometries but required modified Lennard-Jones potentials to accurately reproduce the QM interaction energies. The modified parameter set allows accurate modeling, including free energies, of cation-π indole-choline and indole-trimethylated lysines interactions relevant for protein-ligand, protein-membrane, and protein-protein interfaces.

Search and Subvert: Minimalist Bacterial Phosphatidylinositol-Specific Phospholipase C Enzymes.

Phosphatidylinositol-specific phospholipase C (PI-PLC) enzymes from Gram-positive bacteria are secreted virulence factors that aid in downregulating host immunity. These PI-PLCs are minimalist peripheral membrane enzymes with a distorted (ßα)8 TIM barrel fold offering a conserved and stable scaffold for the conserved catalytic amino acids while membrane recognition is achieved mostly through variable loops. Decades of experimental and computational research on these enzymes have revealed the subtle interplay between molecular mechanisms of catalysis and membrane binding, leading to a semiquantitative model for how they find, bind, and cleave their respective substrates on host cell membranes. Variations in sequence and structure of their membrane binding sites may correlate with how enzymes from different Gram-positive bacteria search for their particular targets on the membrane. Detailed molecular characterization of protein-lipid interactions have been aided by cutting-edge methods ranging from 31P field-cycling NMR relaxometry to monitor protein-induced changes in phospholipid dynamics to molecular dynamics simulations to elucidate the roles of electrostatic and cation-π interactions in lipid binding to single molecule fluorescence measurements of dynamic interactions between PI-PLCs and vesicles. This toolkit is readily applicable to other peripheral membrane proteins including orthologues in Gram-negative bacteria and more recently discovered eukaryotic minimalist PI-PLCs.

A model for hydrophobic protrusions on peripheral membrane proteins.

With remarkable spatial and temporal specificities, peripheral membrane proteins bind to biological membranes. They do this without compromising solubility of the protein, and their binding sites are not easily distinguished. Prototypical peripheral membrane binding sites display a combination of patches of basic and hydrophobic amino acids that are also frequently present on other protein surfaces. The purpose of this contribution is to identify simple but essential components for membrane binding, through structural criteria that distinguish exposed hydrophobes at membrane binding sites from those that are frequently found on any protein surface. We formulate the concepts of protruding hydrophobes and co-insertability and have analysed more than 300 families of proteins that are classified as peripheral membrane binders. We find that this structural motif strongly discriminates the surfaces of membrane-binding and non-binding proteins. Our model constitutes a novel formulation of a structural pattern for membrane recognition and emphasizes the importance of subtle structural properties of hydrophobic membrane binding sites.

Conservation of intrinsic dynamics in proteins-what have computational models taught us?

The intrinsic dynamics of proteins has been suggested to be the most conserved compared to its sequence or structure. As such, the contributing factors to the conservation of dynamics have yet to be determined definitively. Some have suggested that function drives the conservation of protein flexibility, while others have indicated that the overall topology determines protein flexibility patterns. In general, many characteristic features of protein flexibility can be derived from simple coarse-grained models whose success rests on the link between protein local packing density and flexibility. Those models have revealed the evolutionary conservation of protein flexibility and have given us insights into the slow dynamics required for protein function and mechanistic insight into the allosteric effect.

d-Peptides as inhibitors of PR3-membrane interactions.

Proteinase 3 (PR3) is a neutrophil serine protease present in cytoplasmic granules but also expressed at the neutrophil surface where it mediates proinflammatory effects. Studies of the underlying molecular mechanisms have been hampered by the lack of inhibitors of the PR3 membrane anchorage. Indeed while there exist inhibitors of the catalytic activity of PR3, its membrane interfacial binding site (IBS) is distinct from its catalytic site. The IBS has been characterized both by mutagenesis experiments and molecular modeling. Through docking and molecular dynamics simulations we have designed d-peptides targeting the PR3 IBS. We used surface plasmon resonance to evaluate their effect on the binding of PR3 to phospholipid bilayers. Next, we verified their ability of binding to PR3 via fluorescence spectroscopy and isothermal titration calorimetry. The designed peptides did not affect the catalytic activity of PR3. A few peptides bound to PR3 hydrophobic pockets and inhibited PR3 binding to lipids. While the (KFF)3K d-peptide inconveniently showed a significant affinity for the lipids, another d-peptide (SAKEAFFKLLAS) did not and it inhibited the PR3-membrane binding site with IC50 of about 40µM. Our work puts forward d-peptides as promising inhibitors of peripheral protein-membrane interactions, which remain high-hanging fruits in drug design.